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Abstract 

The addition of noise has a regularizing effect on Hermitian matrices. This effect is 
studied here for H = A + V, where A is the base matrix and V is sampled from the 
GOE or the GUE random matrix ensembles. We bound the mean number of eigenvalues 
of H in an interval, and present tail bounds for the distribution of the Frobenius and 
operator norms of H~^ and for the distribution of the norm of applied to a fixed 
vector. The bounds are uniform in A and exceed the actual suprema by no more than 
multiplicative constants. The probability of multiple eigenvalues in an interval is also 
estimated. 


1 Introduction 

It is often the case that disorder has a regularizing effect on the spectrum of an Hermitian 
matrix. Recall the Wegner estimate [33], which express the regularizing effect of diagonal 
disorder, and which is central in the spectral analysis of random operators. The estimate 
was formulated for matrices of the form 


A + (1.1) 

where A is Hermitian and is diagonal with entries independently sampled from a 

bounded probability density p on M. For such N x N matrices, one has uniformly in A: 

E {eigenvalues of {A + in /}] < ||p||ooiV |/| for any interval / C M , (1.2) 

where ||/o|loo is the essential supremum of p, and |/| the Lebesgue measure of I. The following 
related estimate is also valid; 

P{|(A + U‘^''^p-p >t} < for all j = 1,...,A . (1.3) 
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Presented here are somewhat analogous bounds for matrices of the form 
H = /l«y” + 'J/GOE ^ ^Herm^-j^GUE ^ 

the hrst case concerning a real symmetric base matrix perturbed by a random matrix 
yGOE sajjipiefj from the Gaussian Orthogonal Ensemble, and the second case concerning an 
Hermitian base matrix perturbed by a random matrix sampled from the Gaussian Unitary 
Ensemble. The superscripts, which are displayed here for clarity, will often be omitted. 

The invertibility properties of H = A + V are quantihed in several ways: i) tail bounds 
for the distribution of the norm of H~^ip when is a hxed vector, ii) corresponding bounds 
for the Frobenius and operator norms of H~^, iii) a bound on the expected number of eigen¬ 
values of H in an interval. The bounds are uniform in A and exceed the actual suprema 
by no more than multiplicative constants, as can be seen by considering the case A = 0 (cf. 
Section [7]). 


To state the results precisely we first recall the dehnitions of the invariant ensembles. 
These consist of Hermitian matrices of the form 




Vm ’ 


(1.5) 


where X is an iV x iV matrix with independent standard real Gaussian entries in case of 
GOE, or independent standard complex Gaussian entries in case of GUE, and the asterisk 
indicates Hermitian conjugation. In both cases the probability distribution of V is of density 
proportional to 




with respect to the Lebesgue measure on matrices of the corresponding symmetry: real sym¬ 
metric (GOE, with = 1) or complex Hermitian (GUE, with (3 = 2). The distributions are 
invariant under conjugation by the corresponding class of unitary matrices (cf. P, [IHl l25] . 
where various aspects of the invariant Gaussian ensembles are discussed). 


Throughout we write || 9 ?|| for the Euclidean norm of a vector </?, and for a matrix R write 
||i?||F = \/Tr RR* for the Frobenius (Hilbert-Schmidt) norm and ||i?||op = max<^^o II-Rt^II/IIt’II 
for the operator norm. The following pair of theorems states our main results. 

Theorem 1. 

If either: A is an N x N real symmetric matrix, y) G and V is sampled from GOE, 

or: A is an N X N Hermitian matrix, y) G and V is sampled from GUE, 

then the following bounds apply to the matrix H ■.= A + V, with a constant C < oo which is 
uniform in N, A, and ip: 

1. (Eixed vector) for all t > 1, 

p{l|ii-VII>t\/iv||^ll}<k, (1.6) 
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2. (Frobenius and operator norms) for all t > 1, 

< P{||i/-iF>tiV} < j, (1.7) 

3. (Mean density of states) for any interval / C M, 

P [# { eigenvalues of H in I}] < CN\I\ . (1.8) 

The key for the three statements listed in Theorem [1] is the single-vector tail estimate 
fll.hp . In our approach the two other bounds are concisely derived from it. The main technical 
step in the proof of fll.bp is (12.7p of Lemma 12.11 Estimates of similar nature (concentration 
bounds on quadratic forms) have also played a role in the work of Maltsev and Schlein [TTlfTS] 
and Vershynin |3T]. The lemma is proved below by a Fourier-analytic method. 

Because of the similarity between (11.8p and (II.2p of |33], the former bound may be referred 
to as a Wegner-type estimate (though the similarity of the bounds does not extend to their 
derivations). For GUE perturbations of Hermitian matrices a bound of the form (II.8p on the 
density of states, from which (II.7p can be deduced for that case, was also recently proved by 
Pchelin [26], building on the work of Shcherbina |29| . 

For the next statement, we denote, for Borel sets B cM. and Hermitian matrices H: 

ff{B) = A7(B; H) := number of eigenvalues of H in B. 

Theorem 2. Let H = A + V be as in Theorem\^ Then there is a constant C < oo, uniform 
in A and N, such that for every 1 < k < N and every interval / C M, 

p{V(/)>u<^GTT, (1.9) 

Moreover, for any k-tuple of intervals /i, • • • , C M, 


k 

E - l)+ ■ ■ ■ (A7(4) -k + l)+] < 1[{C\If N). (1.10) 

i=i 

Continuing the comparison with bounds which are known for operators with random 
potential, the case /c = 2 of (II.9p is reminiscent of the Minami bound for matrices A + 
with diagonal disorder, as in (II.ip . for which it was established in [20| that 

P {(H -I- has at least 2 eigenvalues in /} < C (||p|loo |7| iV)^ . (l-H) 

This estimate was instrumental in Minami’s proof of Poisson local eigenvalue statistics for 
random Schrodinger operators in throughout the regime of Anderson localization. Ex¬ 
tensions to /c > 2 were subsequently presented by Bellissard, Hislop, and Stolz [2], Graf and 
Vaghi m, and Combes, Germinet, and Klein [6]. In particular, our derivation of Theorem [2] 
has benehtted from the strategy of |6|. 
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Applications The above bounds are useful for a number of problems in the theory of 
random operators, particularly, pertaining to random band and Wegner-type operators, some 
of which will be discussed in |27]. The estimate (II.6p plays a key role in the proof of 
localization at strong disorder for the Wegner A-orbital model, and some of its variants, with 
conjecturally sharp dependence of the localization threshold on the number of orbitals. The 
bound (II.8p enables density of state estimates for a class of models including the Wegner 
orbital model and Gaussian band matrices. Theorem |2] is used to prove convergence of 
the local eigenvalue statistics to the Poisson process in the regime of localization. In such 
applications the sharp dependence of the above bounds on N and t is of value. 

The bounds discussed here are of relevance also from other perspectives. Effects on the 
spectrum of the addition of a symmetric random matrix has been studied in light of appli¬ 
cations in numerical analysis by Sankar, Spielman and Teng |28) (for the case of Gaussian 
random matrices) and by Vershynin [3l] and Farrell and Vershynin [9] (for more general 
distributions). The addition of GOE/GUE and its infinite volume limit were studied by 
Dyson [7] in the context of stochastic evolution, by Pastur [23] in the framework of the limit¬ 
ing eigenvalue distribution for deformed Wigner ensembles. The regularisation effect in the 
infinite volume limit was considered in [3] in the language of free convolution of Voiculescu 

m 


Relation with previous results In presenting some of the related previous results we 
shall invoke the notion of density of states, and the following notation. For an A x iV random 
matrix H, the normalized average i^(-; H) := [A/'(-; H)] (or just u) is referred to as the 

density of states (DOS) measure. When this measure is absolutely continuous, i.e. of the 
form u{d£) = p{S)dS, its Radon density p{S) is called the density of states function. In 
this notation, the bound (11.8p asserts that the DOS measure H) of H = A + V is abso¬ 
lutely continuous, and its density p{S] H) is bounded by a constant independent of N and A. 

While the results presented here focus on bounds which hold uniformly in the base matrix 
A, related questions have been studied for sequences A^ of deterministic Hermitian matrices 
of increasing size for which the density of state measures converge weakly to a limiting 
measure u^oidS). Pastur [23] has shown that in such situations the perturbed operators Aat-I- 
y^OE/GUE more generally + see below) have densities of states which converge 
weakly to a limit which can be determined from z/qo, and which is absolutely continuous of 
density satisfying poo{£) < (c.f. the monograph [25]). 

There are also several results which rely on the Harish-Ghandra formulae da 13115], and 
thus apply to GUE but not GOE perturbation. For the case that A^ are uniformly norm- 
bounded and the perturbation is GUE it is a by-product of the study of local eigenvalue 
statistics by T. Shcherbina [291130] that the Pastur law also holds in total variation distance. 
Thus one can conclude that 

sup supp(£^; Atv 4-< —-I-o(l), A —)■ cxe, 

\\An\\<K £ TT 

This bound is similar to the GUE case of (II.8p . but it requires the deformation to be bounded. 
For the case of possibly unbounded Hermitian matrix perturbed by the GUE, Pchelin proved 
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that 


sup sup sup p{8] Am + i 

N An e 

i.e. our bound (II.8j) on the mean density of states; his argument builds on |29| . 

The above question was considered also in the more general setting obtained by replacing 
GOE/GUE by Wigner matrices for which the entries above the main diagonal are iid 

though not necessarily Gaussian. Vershynin [31] showed that in such case 

sup p|||(/l^ + E^'®)-lop>tiv| < ^ + 2exp(-iV'=^) (1.12) 

with constants Ck, > 0 depending only on K. Vershynin’s result holds under very mild 
assumptions on the matrix entries; an inspection of the proof shows that if the entries are 
themselves regular (for example, have density bounded by C\/N), the estimate holds without 
the term 2 exp(—We also mention that Nguyen [21] showed that for any K > 0 and 
6 > 0 there exists a > 0 so that 

sup P|||(A^ + > iV“| < iV"' . (1.13) 

Upper bounds on the probability of two close eigenvalues were proved by Nguyen, Tao and 
Vu [22]. 

Recently, universality of local eigenvalue statistics for deformed Wigner ensembles was 
studied by O’Rourke and Vu [23], Knowles and Yin |15l Section 12] and Lee, Schnelli, Stetler, 
and Yau |16| . 

Among the results pertaining to Am = 0, that is concerning the density of state of the 
Wigner matrices without these being used as deformations of a base matrix, we mention 
only a few most relevant to the current discussion. 

One of the forms of the Wigner law asserts that if Vm is sampled from a Wigner ensemble 
of dimension N x N then 

p(£; iT") ^ iy(4-f2)+. iV ^ oo (1,14) 

in the weak sense B|25]. In the special cases of GOE and GUE, this may be strengthened 
to uniform convergence |19j . yielding 

sup p{S, ^ jy ^ oo_ (1-15) 

S TT 

This implies a bound similar to (II.8p but for Ajv = 0. 

Maltsev and Schlein [18] proved that the Wigner law (I1.14h holds in the topology of 
uniform convergence in [—2 + 5, 2 — h] (for an arbitrary 5 > 0) for a class of Wigner matrices 
the entries of which obey certain regularity assumptions. Their results imply that (11.151) 
with the restriction \£\ < 2 — 6 holds for this class of Wigner matrices. The paper [18] builds 
on earlier work by Erdos, Schlein, and Yau [8] and Maltsev and Schlein mi, where it was 
shown that there exists an absolute constant G > 0 for which 

p{||(U^'^-T)-iop>tiv}<p t>l. 
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2 Fixed vector bound 


In general terms, the Wegner bound concerns the inverse of a quantity which fluctuates 
due to the presence of random terms in H. For the bound it suffices to focus on the 
ffuctuations resulting from the randomness in the single term However, for the case 

considered here the contribution of any single diagonal term is too small (by a factor of \/iV) 
for the claimed result. Instead, our proof of the bound (II. 6 p . for a given N x N matrix 
A, and a given vector (^, will focus on the fluctuations in ||i7“^(^|| due to the N random 
variables which determine V(p. We start by reducing the claim to a technical estimate whose 
proof will be given separately, in Section [ 6 l 


The real (GOE) case. As the distribution of V is invariant under orthogonal conjugations, 
and we aim at results which hold uniformly in A, we may assume without loss of generality 
that if = ei, the hrst vector of the standard basis in M^. The matrix V has the form 


V = f ^^0 ^ 

V^V 9 Wj’ 


( 2 . 1 ) 


where go ^ M., g G and W is an (N — 1) x (N — 1) symmetric matrix, go, g and W are 

independent, and go is a standard real Gaussian and 5 ^ is a standard real Gaussian vector 
(i.e., with independent entries having the standard real Gaussian distribution). Thus we 
may write 

A _i_y ^ ^ ( y/^9o + 0 {g + by 

Vn\ 9 + b W + D 

for deterministic a G M, 6 G and D an {N — 1) x {N — 1) symmetric matrix. Set 

Q := {W + D)-y 



Inverting using the Schur-Banachiewicz formulae, we obtain 


Vn 


{A + V)-^e,= 


Therefore 
1 


y/N 


||(A + V)-^ei|| = 


V 2 go + a-{g + byQ{g + b) \ -Q{9 + b) J 

+ WQig + b)\y 


(2.3) 


< 


\V 2 go + a-{g + byQ{g + b)\ 

1 


\\Qi9 + b)\\ 


\V^ 9 o + a ~ (s' + byQ{g + b)\ \V^go + 0. ~ {9 + byQ{g + b)\ 

For any deterministic d, and any t > 0, 


(2.4) 


P 


V^9o + d 


>t)-< 




therefore, hrst conditioning on g and Q, one may conclude that 


P 




\V^ 9 o ~ {9 byQ{g + b)\ 2j y/nt 


(2.5) 
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Combining fl2.4p with fl2.5p one arrives at the key bound 


P<^ ||(kl + i/)-^ei|| > tVN\ < ^ + P 


WQid + ^)ll 


\ I ■\/2go + a — {g + byQ{g + b) \ 


> 


( 2 . 6 ) 


For the second term we have the following estimate, whose proof is deferred to Section |6l 


Lemma 2.1. Let Q be a (non-random) non-zero real symmetric matrix, and let g be a 
standard real Gaussian vector of the same dimension. Then, for any real vector b and any 
real number a. 


P 


\\Qi9 + b)\\ 

K^f + bYQig + b)-a\ ~ J ~ f ’ 


t > 1, 


(2.7) 


for some absolute constant C. 


The estimate (II.6p follows in the GOE case, by combining (12.61) with Lemma 12711 (through 
conditioning on go and Q). 

The complex (GUE) case. Here (12.2p is replaced by 


A ^ ( ^0 + a {g + b)* 

V g + b ^ + D 

where go, a G M, 5 ^, 6 G and W and D are (iV — 1) x [N — 1) Hermitian matrices, go, g 
and W are independent, a, b and D are deterministic, go is a standard real Gaussian and g is 
a standard complex Gaussian vector (i.e., with independent entries having independent real 
and imaginary parts, each of which has the normal distribution with mean 0 and variance 
1/2). Following the same steps as in the GOE case one arrives at 



1 11/ ^ , II ^ yi + \\Q{g + bj^ ^ 1 + \\Q{g + ^)ll 

\/]V \go + 0 ^ — {g + b)*Q{g-\-b)\ \go + a — {g + b)*Q{g-\-b)\ 


(2.9) 


where Q := {W + D) 


To conclude the proof via the arguments used in the GOE case, we rewrite the right-hand 
side of (12.9p in terms of a similar expression involving only real quantities. For this purpose 
we consider with the standard basis (cj)^]^ as a vector space over M with the basis 

(ci, iei, 02, ie2, ■ ■ ■ , e^, ie-N) , 

and denote by Q be the 27V x 27V real symmetric matrix which represents multiplication by 
Q in this basis. For a vector v G C^, denote by h G its image under this identification. 
Then 

WQig + ^)ll = WQig + ^)ll (2-io) 

and, using that {g -\- b)*Q{g -\- b) is real as Q is Hermitian, that 

{g + b)*Q{g + b) = {g -\- bYQifg -\- b). (2-11) 
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Thus 


_ 1 + WQjg + ^)ll _^ 1 + WQjg + ^)ll _ ^2 12 ) 

\go + a — {g + b)*Q{g + b)\ a — [g + byQ{g + b)\ 

Note that ^ is a real Gaussian vector whose entries are independent with variance 1/2. In 
order to work with standard real Gaussian vectors we rewrite this expression as 

l + IIW + pil _ V2(l + \\^^(V2g + V2b)\\) ^ 

\go + a — {g + bYQig + b)\ \V^go + \/2a — {V^g + y/2by-^{y/2g + \/26)| 

where y/2g is standard Gaussian. Using (I2.12j) and (I2.13p with (12.9j) allows to finish the 
proof in the GUE case with the same argument as in the GOE case. □ 


Remark 2.2. Note that we actually proved the following stronger, conditional version of 
aiM)-- for (p = ei, the estimate U.6{} holds conditionally on the sub-matrix obtained by 
deleting the first row and column of V. For a general p, this translates to the following 
estimate, which will be of use in the seguel: 




{u*Hv 


u,v F p 


>} 



(2.14) 


with a constant C which is uniform in A, N and p. 


3 Probenius norm bound 


To deduce the Frobenius norm estimate (II.7p from (II.6p . we employ the following principle. 
A similar strategy was employed by Sankar, Spielman, and Teng [281 Proof of Theorem 3.3] 

Lemma 3.1. Let Q be an N x N real symmetric matrix and p be a random vector uniformly 
distributed on the sphere = {^jJ ^ HV’II = !}■ Then 


Proof. By the Ghebyshev inequality, for any real 


p {iiQ</?ii < (-^ii<5<^ir) • (3-1) 

A uniformly distributed vector on can be generated by letting p = ( 7 /|| 5 '|| with g a 

standard real Gaussian vector, for which and 


II 51 II are independent. Thus, 


E [exp (-^IIQt^II^)] = E [exp )] 

^ -E [exp (-^IIQ^IIVII^IP) l||dk< 2 ^] 


r{\\gr<2N} 


(3.2) 


< 2 Eexp 


2N 


\\Qg\? 
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where use was made of the bound P(|| 5 f|p < 2N) > ^ which follows from IE|| 5 '|p = N. 

Let {£j} be the eigenvalues of Q, with which ||Q||| = ^^e distribution of g is 

invariant under orthogonal transformations, and the eigenvectors of Q form an orthonormal 
basis, one gets (using a known Gaussian integral) for any ^ > 0, 


Eexp 


2N' 


WQgf = Eexp 




2N ^ 
i=i 


n ^ ^ 

t=i \/l + jf^j \/l + ;^l|QII 


Juxtaposing the last inequality with fl3.ip and (13.2p . and substituting ^ = 2 t^Q\\^ ~ 

yields 




2exp (^ligil 


1 + 


= 2\/2i 


eexp 


1 -2e^ 


< 5e. 


□ 


We proceed to prove the Frobenius norm estimate (II.7p in the GOE case. Let be 
a random vector distributed uniformly on the sphere and independent of H, and let 

f > 1. Applying Lemma [3.II with Q = and e = we get 


P {11-^ ^IIf > — P [l||ir-i||F>iAf] < 2 E 

t^/N 


^l|F^'"hlF>iAfP { \\H 7^11 — 


tVN 


10 


H 


< 2 E 


P<II^-VII> 


10 


H 


< 2 P<^ \\H-^ip\\ > 


tVN' 

~w 


(3.3) 


Applying now the hxed vector bound (II.Op conditionally on (p to the probability in the last 
term one gets 

20 C 

P{||W-iF>fiV} < (3.4) 

i.e., (II.7p holds in the real (GOE) case. 


A similar argument may be used to establish (11.71) in the GUE case using the following 
complex analog to Lemma 1301 If Q is an W x W Hermitian matrix and is a random vector 
uniformly distributed on the complex sphere, 

S?-'= {^ 6 C": IIV.il = 1}, (3.5) 

then, for all e > 0, 

p|||Qv>ll<^IIQIlFj<5e, (3.6) 

The inequality follows from Lemma 13.11 applied with 2N in place of N by identifying the 
space with as in the proof of the GUE case of (II.Op . This identihcation multiplies 
the Frobenius norm by \/2- 
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4 Bound on the density of states 


We now turn to the density of states bound (II.8j) . 

Let H be the random matrix of Theorem [H Observe that almost surely H has only 
simple eigenvalues, e.g., as its distribution is absolutely continuous with respect to that of 
the underlying invariant Gaussian ensemble (GOE or GUE) and these are well known to 
have this property |19| . 

For a hnite interval /, let be a nested sequence of partitions of I into subintervals 

whose maximal length tends to zero as M —)■ oo. Using the simplicity of the spectrum, 
almost surely: 

y \{H has an eigenvalue in Ij^m O /} ^ 7 ( 1 ^{eigenvalues of iL in /} . (4-1) 

M^oo 

J 

Taking the expectation value and applying the monotone convergence theorem gives 

E [ 7 (/:{eigenvalues of if in /}] = lim > ¥{H has an eigenvalue in Ij^m H /} . (4.2) 

j 


The probabilities on the right may be estimated through the norm bound fll.7p . which implies 
that for any interval J = [£ — e,£ + e] 


F {H has an eigenvalue in J} = P <| || (TT — ||op > yjj f < CN\J\/2 . (4.3) 


Upon summation this yields the claimed density of states bound fll.Sp . 


5 Minami-type bound 


The proof of Theorem |2] proceeds by induction on /c, using an idea of Gombes, Germinet, 
and Klein p. The case k = 1 is exactly the Wegner-type estimate fll.Sp . Thus we assume 
that fll.lOp is valid for a certain k and prove that it is also valid for k + 1. 

In the proof we use the inequality (12.141) . which we restate for convenience. Letting 
H = A+V be as in TheoremlH ip G M^\{0}, denote by the matrix obtained by restricting 
H to the subspace orthogonal to ip, i.e. the {N — 1) x {N — 1) matrix = P^±HP*j^, where 
P^± is the orthogonal projection onto the orthogonal complement of ip. Then (I2.14p asserts 
that 


P 








c 


( 6 . 1 ) 


with a constant C which is uniform in A, N and ip. 

Let H be as in Theorem [1] and hx Ji to be a hnite interval. Let be a random vector, 
independent of H, which is uniformly distributed on the unit sphere in the real case 

or uniformly distributed on the complex unit sphere (®®® (13.5p ) in the complex case. 

Lemma [ 3 . II in the real case or its complex version (13.6p in the complex case, imply that for 
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every non-negative random variable X, measnrable with respect to H, and every G M one 
has 


E[X] < 2E 
= 2E 


X-P ||(i7-£)-V||> 


II (^-^) 


-ii 


10^/N 


H 


X ■ 1 




(5.2) 


Now, let be a nested sequence of partitions of Ji into subintervals whose maximal 

length tends to zero as M —)■ cx). As in Section 01 the monotone convergence theorem implies 
that 


E [Afih)iAfih) - 1)+ • • • (Afih+i) - k)+] 

= lim E [lAt(4Mn/i)>i (A/'(/2) - 1)+ • • • {^f{h+l) - k)+] . (5.3) 
j 

We focus on estimating a single summand in the last expression. Let J C A be an interval 
with midpoint £. The event that Af{J) > 1 coincides with \\{H — T)“^||op > ijj- Applying 


E [lAt(j)>i(AA(J2) - 1)+ ■ ■ ■ (AA(4 +i) - k)+] 


< 2E 


L||(r/-£:)-i||op>j^ 


2 (AA(J2)-l)+---(AA(4+i)-fc)+l 


II(/^-£)-VII>E^ 


< 2E 




1 (AA(J2)-l)+---(AA(4+0-fc)+ • (5.4) 


5|J|\/jV 


Let H^p be as above, then the eigenvalues of interlace those of H, therefore X’ilj) — 1 < 
Thus, 


E {I 2 ) — 1)+ • • ■ (A/’(/fc+l) — /c) + 

< E M{I 2 ] Hp,) ■ ■ ■ {J\f (Jfc+i; Hp) — k + 1)+ 


= E 


Afih; (AA(4+i; H^)-k + 1)+ P II (TT - ^)-Vll > 


5\J\^/N 




(5.5) 


< 5^1 J|X ■ E [AA(J2; Hp)--- H^) - k + 1)+] , 


where in the last inequality we have applied the estimate fl5.ip to the matrix H — £. By the 
invariance of the underlying Gaussian ensemble (GOE or GUE), the {N — l)-diniensional 
matrix 


Hp = 


N 


N -I 


-H, 


V 1 


conditioned on (^, has the form treated in Theorem [T] Thus the estimate (ll.lOp . applied 
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using the induction hypothesis to shows that 


E [Afih; (AA(4+i; H^)-k + 1)+] < J] {Co\IjW N{N - 1)) 

(5.6) 

i=2 

Putting together fl5.4p , fl5.5p and fl5.6p shows that 

/c+1 

E -!)+■■■ (AA(4+i) - k)+] < 10C\J\N X J](C'o|/,|iV) . 

i=2 

Taking Cq > IOC, the theorem follows by plugging the last estimate back into fl5.3p and 
performing the summation. 


6 Ratio of quadratic forms 


Let us recall from Section [2] that the above results hinge on the estimate stated in Lemma [2Tl 
The statement to be proved is that for any (non-random) non-zero real symmetric matrix 
Q, real vector 5, real number a and t > 1, 


P 


WQig + m 


\{g -h byQ{g + b) -a\ 


>t)-< 


C 

T’ 


( 6 . 1 ) 


where 5 ^ is a standard real Gaussian vector and C is an absolute constant. 


That such a bound may hold may be surmised from the observation that 

niQig + bW < C Var[(r? + bYQig + b) - a] 

(uniformly in Q, b, and a), which implies that the denominator of the ratio in fl 6 .ll) fluctuates 
on a scale which is not smaller than the typical size of the numerator. However, more careful 
analysis is needed to take into account the dependence of the two terms and the possibility 
that the denominator has unbounded probability density at small values. 

We turn to the proof of Lemma [2Tl starting with two preliminary claims. The hrst covers 
its rank one case. 


Claim 6.1. If h is a standard Gaussian variable, a, & G M, then 


P 


\h + b\ /8 1 

\{h + by-a\- J - Vvrt’ 


t > 1. 
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Proof. The event > t coincides with 

> \{h + bf - a\. (6.2) 

If a < 0, the probability of this event will only increase if we replace a with 0, thus we 
suppose that a > 0. Then 

\{h + by — a\ = ||h + 6| — \/a[ ■ \ \h + b\ + \/a| > ||h + 6| — ^/a\ ■ \h + b\, 

whence 

The next claim will be used in deriving probability bounds on ratios through estimates on 
the Fourier transform of the joint probability distribution of the numerator and denominator 
(also known as the joint characteristic function). 

Claim 6.2. Let X > 0, Y be a pair of random variables, and 

X(^,h) :=Eexp(i(^X + r/F)). (6.3) 

Then, for any e > 0 and a G M, 

P i > e~^ I < -—- liminf [ dri 

I |F - a| - j - An s^+o J ' 

Proof. The right-hand side of fl6.4p does not change if we replace Y with Y — a, therefore 
we can assume that a = 0. Set 

lx>o, hs{x, y) = h{x, y) exp(-(5a;)a:^ 

and note that 

h{x,y) > 

lai — 

Therefore by the Chebyshev inequality and the Fatou lemma, 

p| ^ I < e^m{X,Y) < eUimMEhs{X,Y). 

Mr I I (5^.+o 

The function hs is continuous and integrable, and its Fourier transform hs is also integrable, 
as follows from the explicit computation below. Therefore, by a version of the Plancherel 
theorem for the Fourier-Stieltjes transform [TTl §VI.2], 

Eh5(X,F) = jj d^dr]hs{^,r])x{^,v), 


h{x,y) = exp 


4e^a; 


d^ 




2 


1+5 


(6,4) 
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where 

hsU, = jj ^< 5 ( 3 ;, y) exp{-i{^x + 7]y))dxdy. 

To compute hs we first fix x > 0 and integrate over y (using a standard Gaussian integral) 


h{x,y) exY>{-ir]y)dy = 


POO 


/ exp 

/ —oo 

[ Ae^x 


dy = 2 \pKxe exp(— 77 ^ 6 ^ 0 ;). 


Multiplying by e ^^x^ and integrating over x, 


hsi^,v) = [ x^ 2 ^^ exp{—x{r]‘^€‘^ + + 5))dx = ^ ^ 


(^ 2 g 2 _|_ 7 ^ _j_ 2 


4+^ 


This implies 

(^) /[ d^dr]hs{^,r])x{^,v) = ^2^X2^^" 


d^dr]{r]^e‘^+ i^ + 6 ) 2 ^x(^,r 7 )- 


Applying the Fubini theorem and taking absolute value, we hnally obtain: 


P 


|h"l 


>e-' 


1 


< -lim inf 

dvr (S^’+o 


dr] 



x{^,v) 

{rfe^ + 7^ + 


□ 


Proof of Lemma \2.1[ Using the symmetry which is built into the assumptions, it suffices to 
establish the bound for diagonal matrices Q = diag(£^i, 82 ,'' ')■ is to prove that 


P 




(6.5) 


where the sums may be restricted to Sj 7 ^ Q (and the probability average is over the inde¬ 
pendent standard Gaussian variables gj). 


We reorder the eigenvalues {Sj) so that 

S^Al + bl)>S^,{l + bl)>£i{l + bl)>---. 

Denote 

fo, ff(l + 6f)<l3E,>.f|(l + '>J) 

r:=ll, £f(l + i?)>i;E,>.£7(l+'>?). £|(l + i>i)<13Ey>2^/(l + '>?) (6-6) 

I 2 , otherwise 


X :='^Sj{gj+ bj)^, Y :='^Sj{gj+ bj)^, x(^,h) := Eexp(7(^A: + r^U)) , (6.7) 

j>r j>l 


14 














where, according to the number of non-zero eigenvalues, X is either identically zero or almost 
surely positive. Observe that 


+ hr- s E \h\\9i+h\ + 'hx. 

i>i i=i 

For the terms in the hrst sum in the right-hand side of (16.81) . Claim IbTTl yields 


( 6 . 8 ) 


P 


\£j\\gj + 6 


|r-a, 

Thus, to prove fl6.5p it suffices to show that 


t>l. 


P 


c 


\Y-a\ 


>t} < —, t>l . 


(6.9) 


( 6 . 10 ) 


If X is identically zero the inequality is trivial. Thus we assume that X is not identically 
zero and note that this assumption entails that Si, 82,^3 7 ^ 0. We now use Claim Wf2\ which 
reduces the task of proving fl6.10p to showing that 


lim inf / dri 
(5'^+0 / 


d^- 




+ 5') 2 +^' 

Noting that a standard Gaussian random variable h satishes 


< C. 


Eexp(ia!(h -1- /3)^) = 


\/l — 2 ia 


= exp 


za 


1 — 2 ia 


, 


we have 


X{^,V) = E 


n 


exp I i ( + bjf + + vSj){gj + hj 

j=l j>r 

irjSj 


exp ( 6T 


.^1 - 2i'r]Sj V 1 - 2,igS. 


n 


: exp ( 6^ 


ij^S] + gSj) 

1 - M^Sf + r]Sj 


j>r -y/l 2i(^Sj gSj) 

For real g, the function xihV) has an analytic continuation to the domain 


2 ma.Xj^rSj ' ’ 


this continuation is given by 


x{^-iS,v) = n 


X 


L1 v^l - 2igSj 
1 


igSj 
^ 1 — 2igSj 


exp ( 6? 


exp (62 


^(1 - 26S]) - 2 iQ " V ' (1 - 26£’|) - 2<: 


( 6 . 11 ) 


( 6 . 12 ) 


(6.13) 


(6.14) 


(6.15) 
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where we set 


(6.16) 


(j ^Sj + r]Sj. 

Due to the assumption that there are at least 3 non-zero eigenvalues, we have: 

d^dr]—^< oo tor 0 < o < 


(e + <52)s 


2 maxjyr Sj 


Thus we may change the contour of integration and apply the dominated convergence theo¬ 
rem to obtain that 


liminf / dri 
< 5'^+0 / 


d^ 




{jfe^ -|- -|- 6') 2 




< dr] d^ 


\x{^-i5,r])\ 
(^2 + 52)! 


(6.17) 


We proceed to prove fl6.1ip by bounding the right-hand side of fl6.17p for a suitable 6. 


Let 


and observe that 


Then from (I6.15p 


z/2 := 




j>r 


s < 


10 maxjyr £] 


(6.18) 


\x{^-iS,r])\ = JJ 


=1 (1 + 4 ? 72 £: 2)4 


- exp ( —2b 


1,2 'I ^3 


X 


n 


1 


^ 1 + A7]‘^8j 

J£]{l-25£^)-2Q 


,>.((1-2<5£:|)^ + 4C|)^ 

Note that, for our choice (I6.18p of v and 5, 

1 1 


- exp b 


^ {1 - 26£f)^ + 4Cf 


< 


< 


exp{l0(5£:2} 


((1 - 26 £fy + 4 : 0 ) " (1 - 25^2)2(1 + 4p2) - ^ ^ 4^2 


and 


6£‘f{l-26£f) \ ( 562^2 


“p I " yi-W+c 


< exp 


3 3 


1 - 2d£] 


< exp {I06b]£f) 


consequently. 


\x{^-i6,r])\ < e JJ 


. 1 (1 + 4r/2T2)4 


exp 


2r]^b]£f 


1 + 4772^2 


n 


3 2 j>r (7 + 4C|)'‘ 


exp 


2bX 


3^3 


1 + 4C 
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Combining this bound with Holder’s inequality yields 


J dr] 

r 

<en 

i=i 


\x{^-iS,r])\ 

di 


dr] 


X 


(^2 + 52)1 (i + 4r72£:2)3 
d^dr] T-r 1 


exp 


Qr]%']£] 


1 + 4772^5 


n 


-^— exp 

(^2 + ^2)J jj. (j 4rt)io=;i 


66? C 


3-r (6.19) 

3 


3 - r 1 + 4C/ 


j=i 


The hrst r integrals satisfy 




1 

f 

f c-u ( 

Qb’jT]'^ \ 

IVdJ 


(1 + 47^2)! 

1 + 47^2 y 


C, 


m^ + \br\)VS 


<C2, (6.20) 


for absolute constants Ci,C' 2 , where the last inequality uses the choice fl 6 . 6 p of r and the 
dehnition fl6.18p of u and 6. 

It remains to estimate I'. An additional application of Holder’s inequality with exponents 


aj = 


sm + m 




shows that 


/'<n 


j>r 


d^dt] 


(^2 + 52)4 


y (1 + 4C|) exp 


26? 


1 + 4C ■ aj 


/-2 

■’i 


( 6 . 21 ) 


j>r 


We proceed to show that each of the Ij is bounded by an absolute constant. Recalling the 
dehnition fl6.16p of Q and changing variables, 


1 

f 

f dCj 

C| 2'>?1 

\Sj\Vsj 

(i+a^^ 

3 

(1 + iCj) 

1 + 4((y oij 


Cs 

f dQ 

1 — 

CM 1 

1 _ 

\£j\VsJ 

(1 + 4cy) 

1 + 4:(j aj _ 


for an absolute constant C 3 > 0. By the choice fl 6 . 6 p of r and fl6.18p of v and <5, 

^ ^^ for all i > r 

4(3-r)a,- - r)£]il + b]) - A ^ 


( 6 . 22 ) 
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whence, splitting the domain of integration into \(j\ < 1 and |Cj| > 1 , 


f dCj 

Cl 26|1 

^ G4-^/o7 ^ C^\8j\ 

/ 3 

/ _|_ 4^2 j 4(3-r)a,j 

1 + 4:(j aj _ 

~ max(l, \hj\) ~ V 


for absolnte constants C 4 , Cs > 0. Plngging the resnlt into (\6.22h and then into fl6.2ip shows 
that /' is bonnded by an absolnte constant. Combining with the bonnds (I6.20j) and plngging 
into fl6.19j) and (I6.17jl . we conclnde that flb.llj) holds, and therefore so does Lemma [2.11 □ 

7 Discussion 

Sharpness of the estimates The key step in onr discussion of the invertibility properties 
of A + V, for a hxed Hermitian, real or complex, matrix A, and a random perturbation V 
sampled from a corresponding Gaussian random matrix ensemble, was the hxed vector bound 
fll.bp . It may be of interest to note that up to multiplicative constant fll.bp is saturated in 
two very different situations: 

1. A = ff (or slightly more generally A = with \£\ < 2). In this case, 

is typically of the order of the contribution of the closest eigenfunction, and for that, 
typically: 

dist(0, spec(iL)) X 1/iV and |((/?, Ti)| x l/\/]V, (7.1) 

where Ti is the eigenfunction of eigenvalue closest to £. 

2 . A = P^, with the orthogonal projection on the space spanned by and 

its orthogonal complement. Perturbation theory allows to conclude that in this case, 
typically: 

dist(0,spec(i7)) X l/\/]V and |((/?, \hi)| x 1. (7.2) 

In both cases ||P“^(p|| is (typically) of the order of the most singular contribution, which is 
|((p, Ti)| dist“^( 0 , specP), and hence 


||P-V|1 x\/]V (7.3) 

up to a random factor whose distribution has \/t tails. However the composition of this 
bound is quite different in the above two cases. 

Note that, while in the above two cases is of the same order, the same cannot 

be said for the density of states at energy 0: it scales as N in the hrst case (i.e. up to a 
constant as fll.Sp L but only as in the second case. 

The Minami-type bound (11.91) is not expected to be sharp since one expects the eigenvalue 
repulsion to result in a higher power on the right-hand side of ( 11 .9p when k >2 (namely, 
in the GUE case and k{k + l)/2 in the GOE case). 
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Weak disorder limit To probe the effects of weak disorder one may consider operators of 
the form; 

Hx,n = An + (7.4) 

with A > 0 a parameter which allows to tune the strength of the disorder. The bounds 
derived here share the property of the random-potential Wegner estimate, that at weak dis¬ 
order the constants degrade at the rate A“^. 

Question; Can the density of states bound for Hx,n be improved in case the base operator 
Ho^n = is itself asymptotically of a bounded density of states? 

(The question is open and of interest also in the original Wegner case.) 

Wigner matrices It is natural to consider extensions of the bounds in Theorem [1] to 
deformed Wigner matrices, about which much has recently been learned [TSlIlG] . The bounds 
cannot hold for any distribution of the entries; in case 1/ is a Wigner matrix with Bernoulli 
entries (uniformly sampled from x/NA = eiel + M 

II (yl -|- 17 )“^||qp —)■ cx) as M —)■ oo, on the event that Vii = 

In particular, for the supremum over N x N real symmetric matrices we have ; 

supP(||(A-M/)"^||op > t) > ^ for any t, 

in contrast to (II.7p . Still, it seems reasonable to expect that bounds analogous to those 
presented in Theorem [T] should hold when the entries of the Wigner matrix are sufficiently 
regular, e.g. with probability densities bounded by \/]V. 
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